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Abstract. We demonstrate the temporal Talbot effect for trapped matter waves 
using ultracold atoms in an optical lattice. We investigate the phase evolution of an 
array of essentially non-interacting matter waves and observe matter- wave collapse and 
revival in the form of a Talbot interference pattern. By using long expansion times, we 
image momentum space with sub-recoil resolution, allowing us to observe fractional 
Talbot fringes up to 10*^ order. 
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1. Introduction 

Interference of matter waves is one of the basic ingredients of modern quantum physics. 
It has proven to be a very rich phenomenon and has found many apphcations in 
fundamental physics as weU as in metrology p?) since the first electron diffraction 
experiments by Davisson and Germer [2j. Matter- wave optics has now developed into 
a thriving subfield of quantum physics. Many key experiments from classical optics 
have found their counterpart with matter waves, for example the realization of Young's 
double slit experiment with electrons [3j, the implementation of a Mach-Zehnder-type 
interferometer with neutrons [4j, or, more recently, the observation of Poisson's spot with 
molecules [5j. The creation of Bose-Einstein condensates (BEG) in 1995 ^ Cj opened 
the door to many more exciting experiments with matter waves, to a large extent in the 
same way as the laser did in the case of classical light waves. 

One remarkable phenomenon in classical optics is the Talbot eflFect, the self-imaging 
of a periodic structure in near field diflFraction. The eflFect was first observed by Talbot 
in 1836 [8j and was later explained in the context of wave optics by Rayleigh in 1881 
[^. When light with a wavelength A illuminates a material grating with period rf, the 
intensity pattern of the light passing through the grating reproduces the structure of the 
grating at distances behind the grating equal to odd multiples of the so-called Talbot 
length I/Taibot = rf^/A. At even multiples of the Talbot length the intensity pattern again 
reproduces the structure of the grating, but shifted laterally in space by half of the 
grating period. In between these recurrences, at rational fractions n/m of Lxaibot (with 
n, m coprime), patterns with smaller period d/m are formed. This effect is known as the 
fractional Talbot eflFect. A necessary requirement for the appearance of the Talbot eflFect 
and its fractional variation is the validity of the paraxial approximation [lOj. Grucial to 
the Talbot effect is the fact that the accumulated phase differences of the propagating 
waves behind the grating show a quadratic dependence on lateral distance or grating 
slit index. 

The first observations of the atomic matter- wave Talbot effect [HI [12j were based on 
setups comprising an atomic beam and two material gratings, where the second grating 
acted as a mask used for detection purposes. The demonstration of the fractional 
Talbot effect with atomic matter waves used the fact that the interference fringes could 
be recorded directly by using a spatially resolving detector [13j. The Talbot effect can 
also be demonstrated with spatially incoherent wave sources by using an additional first 
grating to create spatial coherence according to Lau \\A^ . In this way, an interferometer 
is formed that is made of two or even three gratings. Such Talbot-Lau interferometers 
[15] are now an important tool in atomic and molecular interferometry [l] [161 llZl • 
the context of macroscopic matter waves, i.e. atomic BEGs, the Talbot effect has been 
observed in the time domain by using pulsed phase gratings formed by standing laser 
waves [18j. During expansion after release from the trap the BEG was exposed to two 
short grating pulses separated by a variable time delay and the momentum distribution 
was measured. At a specific delay, this distribution was observed to rephase to the initial 
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one. In essence, the quadratic dispersion relation of freely propagating, non-interacting 
matter waves resulted in a quadratic phase evolution for the diffracted momentum states 
and hence to a temporal version of the Talbot effect. Intriguingly, the Talbot effect is also 
present for interacting matter waves, as we could show in our previous work [19j. The 
momentum distribution of a trapped array of decoupled two-dimensional BECs proved 
to exhibit a regular, time- varying interference pattern. In this case, the quadratic phase 
evolution was driven by the local mean-field interaction that had a quadratic spatial 
dependence reflecting the parabolic shape of the initial density distribution. 

In the present work, we report on the demonstration of the temporal Talbot effect 
using trapped, non-interacting matter waves. Here, the Talbot effect is not driven by 
interactions but by the (weak) external harmonic dipole-trap confinement, leading to 
a characteristic quadratic phase evolution. For our measurements we use as before an 
array of pancake-shaped, two-dimensional BECs in a one-dimensional optical lattice [19j. 
The optical lattice takes on the role of the grating. Cancelling the effect of interactions 
in the vicinity of a Feshbach resonance and decoupling the individual BECs by means of 
a gravitational tilt initiates long-lived Bloch oscillations (BO) in momentum space [20j. 
These are quickly superimposed by a Talbot-type interference pattern in the presence 
of the external confinement. The pattern can be directly connected to the (fractional) 
Talbot effect. In particular, after specific hold times that are multiples of the Talbot 
time, the time-analogue to the Talbot length, a rephasing of the momentum distribution 
can be observed. 

2. Preparation of the initial sample 

We first produce an essentially pure BEC of Cs atoms (no detectable non-condensed 
fraction) by largely following the procedure detailed in Ref. [2T[ [22] . The atoms are 
in the lowest hyperfine sublevel F = 3, tuf = 3 trapped in a crossed optical dipole trap 
and initially levitated against gravity by a magnetic gradient field. As usual, F is the 
atomic angular momentum quantum number, and mp its projection on the magnetic 
field axis. For the present experiments, the atom number is set to typically 6x10^ atoms. 
The trap frequencies in the crossed dipole trap are chosen to be cj^; = 27r x 21.7(3) Hz, 
(jjy = 27T X 26.7(3) Hz, and (jJz = 27t x 26.9(3) Hz. The confinement along the vertical 
axis (z) and the two horizontal axes (x, y) is controlled by two horizontally propagating 
dipole trap beams with beam waists of 46 /xm and 144 /xm and one vertically propagating 
dipole trap beam with a beam waist of 123 /xm. The atomic scattering length and 
therefore the strength of interactions in the BEC can be tuned via a magnetic offset 
field S in a range between as = Oao and as = 1000ao by setting B to values between 
approximately 17 and 46 G using a magnetically induced Feshbach resonance ^23j as 
illustrated in figure[l](a). Here, ao is Bohr's radius. For the initial preparation of the 
sample, we set to positive values, typically between 100 ao and 210 ag. Later, 
is set to zero as discussed below. We gently load the condensed atomic sample into 
a vertical standing wave as illustrated in figure[l](b) by exponentially ramping up the 
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power in the standing wave over the course of about 1000 ms. The standing wave is 
generated by a retro-reflected laser beam at a wavelength of A = 1064.48(5) nm with 
a 1/e- waist of about 350 /xm. We are able to achieve well depths of up to 40£'r, 
where £'r = /l^/c^/(2m) =/i^/(2mA^) = /cb x 64nK is the atomic photon recoil energy. 
Here A; = 27r/A, m denotes the mass of the Cs atom, h is Planck's constant, and /cb 
is Boltzmann's constant. The lattice light as well as the light for the dipole trap beams 
is derived from a single-frequency, narrow-band, highly-stable Nd:YAG laser that seeds 
a home-built flbre amplifler [24j. The maximum output power is up to 20 W without 
spectral degradation. The powers in all light beams are controlled by acousto-optical 
intensity modulators and intensity stabilization servos. 

3. Phase evolution and the Talbot effect 

Our system, the BEG loaded into a ID optical lattice with spacing rf = A/2, can be 
modelled by a discrete nonlinear equation (DNLE) in one dimension [25j , as discussed in 
our earlier work [19]. In brief, this equation can be obtained by expanding the condensate 
wave function from the Gross-Pitaevskii equation, in a basis of wave functions 
centred at individual lattice sites with index j, ^^(z, r^, t) = J^j '^^)- 
Here, z is the coordinate along the (vertical) lattice direction, rj_ is the transverse 
coordinate, and Cj{t) are time-dependent complex amplitudes. The atoms are restricted 
to move in the lowest Bloch band and we can write z) = WQ\z)^±{pj^r±)^ where 

Wq\z) are the lowest-band Wannier functions localized at the j-th site and ^±{pj^r±) 
is a radial wave function depending on the peak density pj at each site [25] . By inserting 
this form into the Gross-Pitaevskii equation and integrating out the radial direction, the 
DNLE is obtained, 
Oc ■ 

= ^(9-1 + 9+i) + ^r('^i>i + ^i'^i- (1) 

Here, J/h is the tunnelling rate between neighbouring lattice sites, Vj =dd j + Vj^^"^ 
describes the combination of a linear potential with force ^ and an external, possibly 
time- varying trapping potential 1/*^^^, and Ej^^{cj) is the nonlinear term due to 
interactions. 

We flrst load the BEG into the vertical lattice and then allow the gravitational force 
to tilt the lattice potential. We thus enter the limit ^ J, in which tunnelling between 
sites is inhibited and the on-site occupation numbers |cjp are constant, determined by 
the initial density distribution. The time evolution of the system is then given by the 
time-dependent phases of all Cj, and the ID wave function ^{q^t) in quasi-momentum 
space q acquires a particularly simple form [26j : 

J j 

= ^^^^^ e-^(^+t)^'^ ^-i{Mj-df-a,^tU-sr)t/h 
j 
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Here, we have assumed that our external potential is harmonic, given by = /5tr(i — 
5)^, where (3t^ = muoldP / 2 characterizes the strength of the potential with trapping 
frequency uj^ along z for a particle with mass m. The parameter 5 in the interval 
[—1/2, 1/2] describes a possible offset of the potential centre with respect to the nearest 
lattice well minimum along the z-direction. For the interaction term aint, the spatial 
dependence is also parabolic, reflecting the fact that we initially load a (parabolically 
shaped) BEC in the Thomas- Fermi regime. In our experiments, the offset 5 is not well 
controlled. It is nearly constant on the timescale of a single experimental run (duration 
of up to 20 s), but its value changes over the course of minutes as the positions of the 
horizontally propagating laser beams generating the trapping potential and the position 
of the retro-reflecting mirror generating the vertical standing wave drift due to changes 
of the ambient conditions. 

The phase evolution in equation ^ has a simple interpretation. The term in 
the exponent linear in j results in Bloch oscillations [271 EHl ED] with a Bloch period 
^Bioch = 2Tih/ {^d). In flgure[l](c) a full cycle of one BO, corresponding to a Bloch phase 
from to 27r, is shown. When restricting ourselves to times that are integer multiples of 
^Bioch this term can be omitted. The nonlinear exponents proportional to lead to a 
dephasing between lattice sites, resulting in a time- varying interference pattern for the 
quasimomentum distribution [19j. In our experiments we have full control over these 
nonlinear terms, not only over /3tr via the external trapping potential, but also over 
the interaction term characterized by a-^^t via the scattering length a^. Our previous 
work [T9j has focused on the role of interactions, whereas in this work we focus on the 
(nonlinear) term caused by the external potential. For this we tune in such a way 
that the term with a-^^t is minimized. Now the phase evolution depends only on the term 
with /3tr. The offset 5 slightly modifles the Bloch period, resulting in a global shift of the 
interference pattern in quasimomentum space when imaged at integer multiples of the 
original TBioch- However, as it is irrelevant for the Talbot effect, we set 5 to zero here. 
By including the simpliflcations and introducing the Talbot time Txaibot = ^/(^^z^^^)? 
equation ^ reduces to 

*(g, t)=Y^ c,(0, q) e--^"^/(2^Taibot) (3) 

3 

with Cj(0, g) = Cj(0) exp(— igjrf). Now the Talbot effect is evident. For times that 
are even multiples of Txaibot the original wave function is recovered, whereas for odd 
multiples the original wave function appears with a shift of hk in quasimomentum space. 
This realisation of the Talbot effect is nearly ideal, since no paraxial approximation is 
needed and since there is no limitation in time due to decreasing wave packet overlap 
fT8]. For fractions n/m of Txaibot, ^ copies of the original wave function with a spacing 
2hk/m appear, corresponding to the fractional Talbot effect. The evolution of the 
quasimomentum distribution as a function of time can be visualized in terms of so- 
called matter- wave quantum carpets [29]. Such a quantum carpet, calculated by solving 



Temporal Talbot effect for trapped matter waves 



6 



equation ([T]) numerically with the parameters typical to our experiment, is shown in 
figure[2] Note that in this case the more simple calculation based on equation ^ leads 
to the same result. However, equation ([T]) gives us more flexibility in relaxing the 
requirements of harmonic conflnement or negligible tunnelling. We plot the distribution 
as a line density plot with white areas indicating high densities. Only times that are 
integer multiples of Tsioch are shown. After a fast spreading of the quasimomentum 
distribution a regular pattern appears at times for which one expects fractional Talbot 
interferences. The number of peaks in the momentum distributions directly represents 
the fraction t/Txaibot- At Txaibot a refocusing to the initial distribution occurs, shifted 
by hk in quasimomentum space. The evolution is then repeated until at 2TTaibot the 
original wave function is recovered. 



4. Experimental realisation 

For the present experiments we choose a lattice depth of SE^^. For lattice loading the 
interaction strength is set to ag = 100 ao and the external trap frequencies are changed 
adiabatically to populate about 40 lattice sites. After loading, we change cj^ to the 
final value. This change is done sufl&ciently quickly (within 3 ms) to avoid a change 
in the initial distribution due to tunnelling, but sufliciently slowly to avoid motional 
excitations along the z-direction. Then, within 0.1ms, we switch off the levitating 
magnetic field gradient to decouple the individual lattice sites and set the scattering 
length to the value near ag = ao that gives minimal dephasing [20] . Note that the point 
of minimal dephasing does not correspond exactly to ao as residual magnetic dipole- 
dipole interactions have to be taken into account [30]. The shift is calculated to be 
about — 0.7ao. After a variable hold time thoid, which typically corresponds to hundreds 
of Bloch cycles with TBioch = 0.575 ms, we switch the levitation field back on in 0.1ms 
and ramp down the optical lattice and the dipole trap responsible for trapping in the 
vertical direction in 0.3 ms. The ramp is adiabatic with respect to the trap frequency of 
the individual lattice sites, ensuring that the atoms stay in the lowest Bloch band and 
thus mapping quasimomentum onto real momentum [31j. Before taking an absorption 
picture we let the sample expand for 80 ms while it remains levitated and thus map 
momentum to real space. The dipole trap responsible for horizontal trapping is not 
turned off immediately, but instead it is ramped down slowly over the course of 50 ms 
to reduce spreading of the sample in the horizontal direction. At the same time, ag 
is kept at the value that gives minimal interactions to avoid broadening of the sample 
in the vertical direction. From the absorption pictures we calculate the momentum 
width Ap as two times the second moment of the momentum distribution along the 
vertical direction. Note that the presence of the horizontal trap during expansion 
leads to additional broadening in vertical direction. This broadening plus some residual 
incoherent background limits the observable values of Ap. Nevertheless, with our ability 
to image the quasimomentum with high resolution [19] we are able to compare not only 
the momentum width but also the substructure in the momentum distribution to theory. 
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Figure[3] shows the measured momentum distribution of the atom cloud at specific 
hold times thoid that are fractions of the calculated Talbot time Txaibot- For this 
measurement we choose a vertical trap frequency of (jJz = 27r x 22.0(2) Hz, which 
gives Txaibot = 555(10) ms. Figure[3](a) shows the absorption images as density plots 
(white areas indicate regions with high density), while figure[3](b) plots the horizontally 
integrated densities from the corresponding images of figure[3](a). Initially, the 
momentum distribution is singly peaked, as expected for a non-dephased BEC. After 
a rapid coherent dephasing (corresponding to a rapid broadening of the momentum 
distribution, not shown here) regularly structured patterns appear. The number of 
peaks within the first Brillouin zone [—hk^+hk] corresponds exactly to the fraction 
^hoid/Txaibot, ^s cxpcctcd from the theoretical considerations. A small fraction of the 
atoms is detected outside the first Brillouin zone, likely caused by imperfections in the 
mapping of quasimomentum onto real momentum. Finally, at the Talbot time, the 
momentum distribution rephases again to the initial distribution. In general, apart 
from an overall shift of each individual distribution in quasimomentum space due to 
variations in 5 as discussed below, we find very good qualitative agreement with the 
results of the calculation shown in figure[2j 

Figure[4] illustrates the effect of S on the observed patterns in quasimomentum 
space. For two diflFerent hold times thoid = Txaibot and thoid = ^Taibot/S, absorption images 
for several individual experimental realisations and the corresponding horizontally 
integrated densities are shown. The expected single- and double-peaked momentum 
patterns are reproduced from one experimental realisation to the next, but they 
experience a varying shift in quasimomentum space. As a consequence of the periodic 
structure of quasimomentum space, a peak that is located near one edge of the Brillouin 
zone also reappears at the opposite edge. The maximum possible shift of the pattern 
in quasimomentum space due to 6 increases with hold time and is calculated to be 
±hk X thoid/Txaibot. This is why the patterns shown in figure|3| e. g. at thoid = ^xaibot 
or at thoid = rxaibot/2, agree with the calculated patterns only modulo the shift in 
quasimomentum space. Note that, alternatively, we could have chosen to present in 
figure|3] selected patterns from a sufficiently large sample of measurements, e.g. the 
one from experimental run 4 for thoid = ^xaibot or the one from experimental run 1 for 
^hoid = rxaibot/2 shown in figureQ 

A simple quantitative comparison between experiment and calculations can be done 
by considering the time evolution of the momentum width Ap. The distribution of 
this quantity across several experimental realisations is evidently sensitive to the de- 
and rephasing of the matter wave. In fact, we can relax the choice of the Bloch 
phase and allow its value to be random. For example, for a non-dephased BEC the 
momentum width Ap is measured to range from {Ap)^[^ ^ 0.6 hk^ corresponding to 
the singly peaked momentum distribution, to (A]9)max ^ 1.7 hk when the momentum 
distribution is evenly peaked at both edges of the Brillouin zone (e. g. at half the first 
Bloch period, see figure[l](c)). For a completely dephased sample corresponding to a 
uniform distribution over the first Brillouin zone we measure a value of Ap ^ 1.25 hk. 
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Accordingly, the range for the momentum width Ap at a given hold time thoid shows a 
distinct behaviour as a function of thoid, in particular indicating the revival at Txaibot 
by maximizing the difference L^Ap = (A]9)max — (A]9)min between the extrema of Ap. 
Figure[5](a) shows (A]9)max and {Ap)^i^ as a function of thoid as calculated from equation 
Q. Initially and at Txaibot the extrema lie far apart (at these times the calculation gives 
values for {Ap)^[^ that are close to zero in accordance with the fact that the momentum 
width is determined only by the spread in position space, which is large), whereas 
at intermediate times the difference is drastically reduced, only increasing slightly at 
rational fractions of thoid/^Taibot- In figure[5](b) we plot the measured momentum width 
extrema. These are determined from samples of 10 single measurements at each chosen 
value for thoid- The initial rapid collapse agrees well with the fact that the sample 
dephases. Then, near the calculated value for Txaibot, a clear increase in D^^ can be 
seen. The difference recovers almost completely to the initial value. We attribute the 
slight reduction to additional dephasing mechanisms not included in our simple model 
as discussed below. 

The behaviour of L^Ap offers a simple method to test the dependence of the Talbot 
time Txaibot on the vertical trap frequency uj^. Evidently, L^Ap has a maximum at Txaibot- 
Figure[6](a) shows the momentum width Ap in the vicinity of the calculated Txaibot, here 
for a specific trap frequency of cj^ = 27r x 26.9(2) Hz. Again we evaluate 10 experimental 
realisations for each hold time and select (A]9)niax and {Ap)^i^ to calculate i^Ap- We 
locate the position of its maximum by a simple gaussian fit, as shown in figure[6](b). 
We then vary uOz and determine Txaibot accordingly. In Figure[6](c) Txaibot is plotted as a 
function of cj^. The experimental values are in excellent agreement with the calculated 
values for the Talbot time according to Txaibot = h/ (muj'^cP). 

We finally discuss the main limitations for our experiment. We believe that the total 
number of subsequent revivals that we can observe (we detect up to 4 revivals) is mainly 
limited by three-body loss and by the anharmonicity of the trapping potential. Three- 
body loss heats the two-dimensional BECs residing at each lattice site. This leads to a 
loss of phase coherence and thus decreases the visibility of the revivals. Perhaps more 
interestingly, the anharmonicity of the trapping potential along the vertical direction 
causes deviations from the quadratic phase evolution required for the Talbot effect. In 
order to test this effect we generate the vertical trapping potential with a more tightly 
focused dipole trap beam, which enhances the effect of anharmonicity. We then observe 
non-perfect Talbot revivals followed by subrevivals as can be seen in figure |7|(a). This 
is in qualitative agreement with calculations shown in figure[7|(b), for which the real 
gaussian shape of the trapping potential instead of a simple harmonic one has been 
used. The full calculated time evolution of the momentum distribution is shown in 
figurelTFc). The distortion of the matter- wave quantum carpet can clearly be seen. 
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5. Conclusion 

We have demonstrated the temporal Talbot effect with trapped, non-interacting matter 
waves. High resolution imaging in quasimomentum space allows us to resolve Talbot 
fringes up to the 10th order. We have tested the dependence of the Talbot time on the 
strength of the confinement and have found very good agreement with the calculated 
value. We find that the interference pattern is sensitive to the anharmonicity of the 
trapping potential. In principle, the detailed structure of the interference pattern 
and the precise revival times are sensitive probes for force gradients and interactions 
between atoms. The weak magnetic dipole-dipole interaction, for example, has recently 
been investigated in the context of matter- wave interferometry [30j. Matter- wave 
interferometry in the Talbot regime could potentially be used to examine in detail the 
effect of the long-range nature of such an interaction. Similarly, a spatially dependent 
force like the Casimir-Polder force [32[ [33l [34] near a surface could be investigated 
through its infiuence on the Talbot interference pattern. 
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Figure 1. (a) Magnetic-field dependence of the scattering length for Cs atoms 
in F = 3, mi? = 3: Wide tunability is given by a broad magnetic Feshbach resonance 
with a pole near — 11 G (not shown), leading to a region with attractive interaction, a 
zero crossing at about 17 G, and a repulsive region above [23]. Two narrow Feshbach 
resonances can be seen in the vicinity of 50 G. (b) Experimental configuration: A 
vertically-oriented standing laser wave creating a stack of pancake-shaped traps is 
intersected by two horizontal laser beams, (c) Bloch oscillations: Time series in steps 
of about 57 /is showing the quasimomentum distribution over the course of one Bloch 
cycle. 
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Figure 2. Calculated BEC-based temporal Talbot effect, (a) Momentum distribution 
as a function of the hold time thoid starting from the initial BEC to the first revival at 
Txaibot foi" a pure harmonic potential. White areas indicate a high occupation of the 
respective momentum state, (b) Same as in (a) from thoid = ^Talbot to thoid = 2TTaibot- 
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Figure 3. BEC-based temporal Talbot effect - experiment, (a) Series of absorption 
images after 80 ms of expansion, showing fractional Talbot fringes of different order 
in momentum space, starting from the initial momentum distribution of the BEC 
after two BOs (left), followed by the 10*^ order at Txaibot/lO, 9*^ order at Txaibot/Q, 
etc., down to the 0*^ order at the Talbot time (right). Note that the time axis is not 
linear. White areas indicate higher density, (b) Horizontally integrated density profiles 
obtained from the absorption images shown in (a). Note that, e.g. for Txaibot/lO, the 
outermost momentum component appears twice, i.e. at both edges of the Brillouin 
zone. 
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Figure 4. Variations in the momentum distribution between successive experimental 
realisations for long hold times, (a) Absorption images of five individual experimental 
realisations with thoid = ^Talbot- White areas indicate higher density, (b) Horizontally 
integrated density profiles obtained from the absorption images shown in (a), (c) 
Absorption images of five individual experimental realisations with thoid = ^Talbot /2. 
(d) Horizontally integrated density profiles obtained from the absorption images shown 
in (c). Note that, in addition to the random shift in quasimomentum space caused by 
(5, effects of horizontal dynamics, especially fragmentation and density variations along 
the horizontal axis, can be observed. 
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Figure 5. Talbot revival as evidenced by the spread of the momentum width Ap. 
(a) Calculated (Ap)max (blue diamonds) and (Ap)min (black circles) as a function of 
^hoid ill units of Txaibot- (b) Measurement of (Ap)max (blue diamonds) and (Ap)min 
(black circles) as a function of thoid in units of Txaibot for a vertical trap frequency 
of cj^ = 27r X 22.0(2) Hz. The extrema are determined from a sample of 10 single 
experimental realisations for each value of thoid- 
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Figure 6. Talbot time Txaibot as a function of the external confinement strength, 
(a) Momentum width Ap in the vicinity of the expected Txaibot foi" a dipole trap 
frequency of ujz = '27rx 26.9(2) Hz for 10 single experimental realisations (black circles). 
The extrema (Ap)max and (Aj9)min are indicated as red diamonds, (b) Calculated I^Ap 
for the measured extrema in (a). The solid line represents a gaussian fit, from which 
^Talbot is derived, (c) Dependence of Txaibot on the trap frequency Uz =0Jz/{'^'7^)^ 
The black (blue) circles (diamonds) represent measurements for which the external 
harmonic trap is generated by the dipole trap beam with 46 jam (144 jam) beam waist. 
The solid line gives the calculated values for Txaibot- The vertical error bars are the 
1(7 uncertainty of the maximum position of the gaussian fit as shown in (b). The 
horizontal error bars are equal or smaller than the symbol size. 
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Figure 7. Effect of the anharmonic trapping potential on the momentum distribution, 
(a) Momentum width Ap in the vicinity of the expected Txaibot foi" a dipole trap 
frequency of 27r x 31.1(2) Hz for 10 single experimental realisations (black circles). The 
vertical dipole trap is created by the more tightly focused dipole trap beam with a beam 
waist of 46 /im. The extrema (Ap)max and (Ap)min are indicated as red diamonds, (b) 
Calculated (Ap)max and (Ap)min in the vicinity of the expected Txaibot for the same 
experimental parameters as in (a). For the trapping potential the real gaussian shape 
of the dipole trap is used, (c) Full calculation of the momentum distribution as a 
function of the hold time thoid using the same parameters as in (b). 
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